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We study, within the Ginzburg-Landau (GL) theory of phase transitions, how elastic deformations 
in a supersolid lead to local changes in the supersolid transition temperature. The GL theory is 
mapped onto a Schrodinger-type equation with an effective potential that depends on local dilatory 
strain. The effective potential is attractive for local contractions and repulsive for local expansion. 
Different types of elastic deformations are studied. We find that a contraction (expansion) of the 
medium that may be brought about by, e.g., applied stress leads to a higher (lower) transition 
temperature as compared to the unstrained medium. In addition, we investigate edge dislocations 
and illustrate that the local transition temperature may be increased in the immediate vicinity of 
the dislocation core. Our analysis is not limited to supersolidity. Similar strain effects should also 
play a role in superconductors. 

PACS numbers: 



I. INTRODUCTION 

Superfluids flow without resistance. The existence of 
superfluidity raised the possibility of supersolida^L solids 
in which superfluidity can occur without disrupting crys- 
talline order. Long ago, Chester^ theoretically demon- 
strated the possible existence of a supersolid. If super- 
solids exist, a natural contender would be solid Helium. 3 
Recent torsional oscillator experiments^ on solid 4 He 
pointed to supersolid type features and have led to a 
flurry of activity. In the simplest explanation of the 
experiments,^ a portion of the medium becomes, at 
low temperatures, a superfluid that decouples from the 
measurement apparatus. Such a "Non Classical Rota- 
tional Inertia" (NCR!) effect is known to exist in su- 
perfluid liquid HeliurrP^ which was probed with simi- 
lar techniques EEI Experimental results suggest the ab- 
sence of superfluid features in ideal crystals with no grain 
boundaries.^ Currently, it is not clear if an NCRI lies 
at the core of the recent experimental findings in solid 
4 He. For instance, the required condensate fraction ad- 
duced from a simple NCRI-only explanation does not 
simply conform with thermodynamic measurements^. 
Rittner and Reppy- 11 ' discovered that the putative super- 
solid type feature is acutely sensitive to the quench rate 
for solidifying the liquid. Aoki, Keiderling, and Ko- 
jima reported rich hysteresis and memory effects^ sim- 
ilar to those occurring in glasses 21 . The torsional os- 
cillator findings can arise from material characteristics 
alone EMHini j n particular, the thermodynamics and tran- 
sient dynamics of distributed processes in amor phou s or 
general non-equilibrated solids can currently fi"B 1 1 1 3 1 b_ 
served results. Indeed, later numerical results point to- 
wards such a possibility.^ Notably, recent experimen- 
tal results^ agree with an earlier suggested theory con- 
cerning such transient dynamics.^ The presence of non- 
uniformity in 4 He is also suggested by a criterion compar- 



ing the change in dissipation vs. relative period shift in 
torsion oscillator P21 It may well be that theseglassy and 
superfluid effects are present in solid Helium! 2 ^ An inter- 
esting question concerns the coupling between elastic de- 
fects such as dislocations and superfluid type features.^ 
The coupling of the supersolid transition to impurities 
was discussed in Ref. 2$\ The coupling between superflu- 
idity and elasticity in supersolids and how this may lead 
to a strain-dependent critical temperature was discussed 
in Refs. 1261271 The viable existence of supersolid phase 
is not confined to solid Helium. Other contenders for the 
supersolid state include cold atoms in a confining optical 
lattice!^ There has been much work examining superso- 
lidity in spin systems as well, see, for example, Ref. UHl 
Supersolids constitute a fascinating state of matter and 
appear in a host of systems. 

This article focuses on th e cou pling between nanoscale 
structure and supersolidityP^Zl As is well appreciated, 
elastic strain may fundamentally affect local and meso- 
scopic electronic, magnetic and structural properties. 
There is ample evidence for significant coupling amongst 
the electronic degrees of freedom with the latti ce dis tor- 
tions in cuprates, manganites, and ferroelectrics! 3 ™ 31 l The 
central thesis of this work is that elastic distortions may 
alter the supersolid behavior. As we will elaborate later 
on, in, e.g., a cylindrical torsional oscillator geometry in 
which the boundary of the solid is elastically deformed 
so that it undergoes a supersolid transition at a higher 
temperature than the bulk, a fraction of the boundary 
will become a supersolid leading to a partial decoupling 
of the bulk from the torsional oscillator chassis and a 
consequent reduction in the period. 

In this work, we will employ a Ginzburg-Landau (GL) 
theory to study the influence of clastic strain on superso- 
lidity. As we will show, the Eulcr-Lagrange equations for 
the GL free energy result in an effective Schrodinger type 
equation. We find the lattice distortion acts as an effec- 
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tive potential for the supersolid order parameter. Solv- 
ing the resulting effective Schrodinger-type equation, we 
find our main results: (1) a contraction (expansion) of 
the lattice edges leads to an increase (decrease) in the lo- 
cal supersolid transition temperature; (2) elastic defects, 
such as dislocations, lead to similar effects. 

Although our motivation is the analysis of super- 
solids, all of our calculations within the GL frame- 
work are ident ical fo r non-uniform elastically strained 
superconductor^U22l an( j i eac [ to the- same general con- 
clusions, which we will derive in this work. The case 
of uniformly strained superconductors has been investi- 
gated in detail in myriad experiments, starting from Ref. 
l33l and many works since™ It was found in these works 
that uniform hydrostatic pressure can increase the su- 
perconducting transition temperature. The influence of 
pressure on the superconducting temperature has also 
been investigated in numerous theoretical treatments, 
e.g., Refs. 35 3H1 Our GL formulation and Schrodinger- 
type equation give rise to an increase of the supercon- 
ducting transition temperature under applied pressure. 

The outline of the paper is as follows: in Section [TTJ we 
set up the general GL framework for our investigations. 
We illustrate the connection between the Euler-Lagrange 
equation and the Schrodinger equation. In the sections 
thereafter, we focus on particular lattice distortion pro- 
files to determine the change in the local supersolid tran- 



sition temperature. In Section III we examine the influ- 
ence of a boundary edge contraction, and in section |IV 
we study the opposite case of a boundary edge expansion 
In section W\ we analyze the case of an edge dislocation 



We summarize our findings in section VI 



II. GENERAL FRAMEWORK 



We study the GL free energy density 
F(r) = a(T)\ij)\ 2 + ^| 4 + c|V^| 2 + A(r)|V>| 2 , (1) 

where T is the temperature, b and c being positive con- 
stants, tp the (complex) supersolid order parameter, and 
A(r) a position dependent function that captures the cou- 
pling of the order parameter to elastic strain as we elab- 
orate on below. The prefactor b in Eq. ([I]) is positive and 
depends only on the density of the crystal, as well as on 
defect densities.^ For temperatures T < T c , the coeffi- 
cient a(T) is negative enabling a non-zero tp to minimize 
the free energyP^The condition a(T c ) = determines the 
transition temperature T — T c below which supersolidity 
onsets.^ The third term in Eq. relates the free energy 
with the magnitude of the gradient of ip, as in a domain 
wallP^The difference between the free energy of a normal 
crystal and a displaced crystal appears in the last term. 
For a crystal whose constituents i undergo a distortion 
from an ideal unperturbed configuration R to a shifted 
configuration R' due to the application of stresses, we 
set Ui = R 1 ; — Ri and take the continuum limit wherein 



we replace i by the continuous coordinate r. In the up 
and coming, the Greek indices 7, 5 will denote the spa- 
tial components (e.g., 1*7=1,2,3 will denote the Cartesian 
components of the displacement u at site r). In general, 
a linear coupling of the form ajgu^s\tp\ is allowed be- 
tween the linear order strain tensor u^s — ^(dyiis + dsUy) 
(where u is the ela stic di splacement)^ and the supersolid 
order parameter y> P^ 4 ° l In what follows, we will consider, 
for simplicity, the case in which the displacement occurs 
only along one Cartesian direction. Allowing for general 
displacements does not change our conclusions. For uni- 
directional displacements, the coefficient of the last term 
111 Eq. can be expressed as a dilatory strain 



X(r) =dV-u(r), 



(2) 



where d is a positive constant and u(r) is the displace- 
ment field. The sign of d is chosen such that the free 
energy of Eq. ([!]) is lowered on introducing vacancies. 
The vacancy density scales with [— (V ■ u)] (whereas the 
interstitial density scales with [(V-it)]). Ions in the vicin- 
ity of a vacany will have an inward displacement towards 
its location whereas ions in the vicinity of a interstitial 
will be pushed outwards. Eq. ([2| and the free energy are 
functions of the strain tensor and thus symmetric under 
spatial reflections under which r — > — f and u — > —it. In 
bulk linear elasticity, the local strains scale, as in Hooke's 
law, as the pressure divided by the elastic moduli. In the 
following sections, we will consider the strain fields asso- 
ciated with various cases. 

As noted earlier, our GL theory of Eq. ([I]) also de- 
scribes a (singlet) superconductor with an order param- 
eter ip in the presence of elastic strains on which we 
comment below. In a charged crystal (of unit cell vol- 
ume when undeformed), under applied elastic stress the 
electric field couples to the local charge density (which 
deviates from that of the undeformed crystal by an 
amount [—(V • u)] (and whose volume trivially scales as 
(1 + V ■ u))). For superconductors, A(r) — > e*(p(r) w ith <b 
the electrostatic potential and e* an effective chargeP^E] 
With Eq. (|2), the last term in Eq. Q is a general 
isotropic coup ing between the strain and the supersolid 
(superconducting) order parameter. In the cases that 
we will examine the displacement u will occur along one 
Cartesian direction (u will have only one component). 
Furthermore, in the first two cases that we will detail 
below (contraction and expansion along an edge), this 
displacement field will vary only along one Cartesian di- 
rection and will be uniform along all other orthogonal di- 
rections. Consequently, the coupling A will depend only 
on one Cartesian direction: A = A(z). In the last case 
discussed in this work, that of an edge dislocation, the 
displacement field (and consequently the coupling A) will 
depend on two directions. 

To find the ground state of such crystal, we want to 
minimize the free energy. The variational derivative of F 
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with respect to ip* leads to the Euler-Lagrange equation 
SF 



5ip* 



a{T)ip + b\ip\^ip - cV^V + A(r)^ = 0, (3) 



with an identical (complex conjugated) equation for 
SF/Sip = 0. In situations in which a weakly first- or 
second-order supersolid transition occurs, we may, in the 
vicinity of the transition (where ip is small) omit the cu- 
bic term in Eq. ([3|, 42 and the variational equation may 
be recast as 



cV 2 ip 



\{r)ip = -a(T)ip. 



(4) 



Eq. Q is a Schrodinger type equation with c = H 2 /2m 
and a[T) = —E with E the energy and m a mass. Solving 
for the eigenvalue E — —a(T) enables us to extract the 
transition temperature. Generally, a shift in the transi- 
tion temperature results from the coupling to the elastic 
displacements. 

The gradients of u as embodied in A(r), take on the role 
of a potential energy in the effective quantum problem 
for the "wavefunction" ip. We briefly comment that the 
case of uniform pressure corresponds to a constant (V-u) 
and thus to a constant effective potential A(f) = const. 
Applied to the analysis to be presented below to more 
complicated cases, such a uniform shift of the potential 
energy (and thus to the eigenvalues E) leads to a constant 
shift in the value of a(T) at the transition point. As a(T) 
is monotonic in temperature, for d > (d < 0) this leads 
to an increase (decrease) in the transition temperature 
for a uniform contraction (V •u < 0) as it may indeed oc- 
cur under uniform applied pressure in superconductors 
for which for an increase or decrease of th e supercon- 
ducting T c appear for different systems? 33 * 34 ! The case of 
a spatially uniform dilatory stress is a particular simple 
limiting form of the more general non-uniform elastic de- 
formations that we discuss in this work. 

In the remainder of this work, we will examine the 
solutions of Eq. Q for various non-uniform elastic dis- 
placements u. In particular, we will examine the strain 
fields associated with a contraction of the sample bound- 
aries, an expansion of a boundary edge, and the strain 
profile associated with an edge dislocation. 



III. CONTRACTION OF BOUNDARY EDGE 

Consider a crystal with a side of length L along one 
of the Cartesian directions (the coordinate values corre- 
sponding to this side are in the range L/2 > z > —L/2). 
We consider a contraction in which near the two edges, 
the lattice sites are most displaced from their equilib- 
rium positions, see Fig. [Tj Such a contraction may, e.g., 
be brought about by applying stress (along opposite di- 
rections) on the two edges of the system. Alternatively, 
a shock wave or generation of coherent phonon propaga- 
tion by ultrafast pump-probe spectroscopy may be used 



to create local density modulations resulting in nonuni- 
form strain near the edges of the sampled As we will 
show, the displacement at the edges leads to a change 
in the local transition temperature. A displacement field 
describing a contraction along the z direction is given by 

'] for \z\ < L/2 
for \z\ > L/2 ' 
(5) 

The displacement thus occurs in some finite region (of 
scale k) about the edges. Uq the maximum displacement, 
and with no displacement along the x or y directions, 

U X = Uy = 0. 



by 



The corresponding effective potential of Eq. ([2| is given 



2uod 



[z 

V 2' 



(z+|) e -(*+*) 2 /* 2 



( 6 ) 

For \z\ > L/2 (points outside the crystal), the supersolid 
order parameter -0 = and in Eq. Q the effective po- 
tential A = oo. For small deformations, this attractive 
potential leads to the appearance of a weak bound state. 
For z > 0, L/2 > (L/2 - z) > kyj\n(2u d/k) = e/2, 
the effective potential tends to zero, and the bound state 
wavefunction is of the form ip ~ exp[«;(z — L/2)]. A 
similar form is attained near the point z — —L/2. The 
value of k and thus of the bound state energy E = —ck 2 
can be computed in the standard way by integrating the 
Schrodinger equation once in a region of width e across 
the point z = L/2 in an extension of the problem to 
z > L/2 in which the potential is symmetrized about the 
point z = L/2. As \E\ <C |A(,z)| in the narrow region 
near the edges, Eq. Q reads 



- 2k = 



dip 
dz 



L/2+e/2 



1 



£/2-e/2 C JL/2-e 



L/2+e 



\{z)dz. (7) 



Ignoring exponentially small corrections, we attain that 
the bound state energy is 



E= — 



4c 



(8) 



We will now employ the value of E to determine a change 
in the transition temperature. Within the GL theory, 
a(T) ~ a(T — T ( ?) near the transition temperature, where 
Tj? is the unaltered transition temperature and a > is 
a constant. Writing a + E = a(T - T c e// ) where T^ ff is 
the effective transition temperature, we have 



with 



AT r 



-?2 ,2 



4ac 



(9) 



(10) 



In other words, the region near the contracted edges has 
a higher transition temperature into the supersolid state 
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distance z 

Figure 1: (Color online.) The displacement corresponding to a contraction near the edges. For clearly illustrating 
the contraction, we depict exaggerated displacements. In this sketch, the displacement (black solid) is given by 
Eq. (J5j) with L = 20, uq = 0.5, d = 1, and k = 1 where the lattice constant set to unity. The effective potential (red 
dashed) is given by Eq. The large (exaggerated) value of the displacement u is chosen to vividly illustrate the 
contraction. 



than the bulk. Generally, the maximal displacement in 
Eqs. ([§ and @ can be of order uq ~ 0.1 lattice con- 
stants as set by the Lindemann criterion of melting in 
most materials (or of Uq ~ 0.2 in solid 4 He and poten- 
tially other quantum solids). 44 In Eqs. and |2]), the 
parameters c,d — 0(1). We estimate a small enhance- 
ment of the transition temperature in the surface region. 
For parameters a = 1/Tj?, d = 1,Uq =0.1, and c — 1, we 



find from Eq. ( 10 1 a small enhancement compared to the 
bulk transition, AT C = 2.5 x 10" 3 T C °. 

The effect of this shifted transition temperature is that, 
when a sample of contracted 4 He is cooled down, the re- 
gion near the edges would turn into supersolid at a higher 
temperature than the bulk of the crystal. Perusing the 
form of the supersolid order parameter ip, and Eq. Q, we 
see that ip drops exponentially away with the boundary 
with a penetration depth I = 2c/[duo). With our previ- 
ous estimates for parameters in Eq. ( |10[ ) combined, we 
find that the penetration depth i ~ 20 lattice constants. 

Returning to the NCFd^HH briefly discussed in the 
introduction, if the entire sample is rotating before the 
transition to a supersolid phase occurs, at some temper- 
ature higher than the normal transition to supersolid of 
bulk helium but low enough to make the edges become 
supersolid, the supersolid component in the edges will 
partially decouple from the bulk rotation. This situation 
is depicted schematically in Fig. pi). 

We now expand on the relation between the local value 
of the supersolid order parameter tp and the local effective 
transition temperature. Eq. Q holds for all locations 
—L/2 < z < L/2 (and trivially, of course, on any local 
segment within this region) . We earlier solved Eq. Q to 
find that the supersolid order parameter decays exponen- 




Figure 2: Under compression of the edges of the 
torsional oscillator, the rim attains a supersolid 
component at a higher temperature than the bulk does. 
On cooling down to this temperature, the supersolid 
fraction of the rim partially decouples from the bulk 
and outer chassis. 



tially with a decay distance I away from the boundary 
points z = ±L/2. Thus, deep within the bulk, the super- 
solid order parameter was zero. We may examine Eq. Q 
locally (with a local effective potential A(z)) in order to 
see when we may attain a finite supersolid order param- 
eter tp(z) at general locations z away from the boundary. 

As the displacement only occurs near the edges, and 
since the change in transition temperature is the result 
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of the displacement, it is reasonable to assume that the 
change in transition temperature can only be detected 
in the region near the edges. For the region inside the 
crystal (far from the edge), the transition temperature 
should remain unaltered. Based on the observation, the 
transition temperature as a function of the z-axis (the 
axis parallel to the length of the crystal) could be de- 
scribed as 



t: 



eff 



(*)=!? + /(*)( 



4ac 



(11) 



where f(z) is a function that rapidly varies from 1 at the 
boundaries z = ±L/2 to zero for positions removed from 
the boundaries. An example is provided by 



ff z \ = e -(z-L/2) 2 /k 2 + e -(z+L/2f/k 2 



(12) 



A contour plot of is depicted in figure [3J for L = 20 
(we set the lattice constant to be unity), k = 1,T% = 2 

d 2 u 2 

K, u Q — 0.5 and -j^- = 0.2 K. The large (exaggerated) 
value of the displacement u a is chosen to lucidly illustrate 
a contraction as in, e.g., Fig. [T]and its effect. 



IV. EXPANSION OF EDGE BOUNDARIES 

The situation of the expansion near the edge bound- 
aries is schematically shown in Fig. [4] As in the case of 
contraction, this may be physically brought about by ap- 
plying opposite stresses (e.g., shear stresses) on the two 
boundaries of the system. In an annular geometry sim- 
ilar to that in Fig. [2] an expansion may result by a 
difference in pressures between the inner and outer parts 
of the cylinder. A typical displacement field u is, in this 
case, given by 

f Ve-^/ 2 ) 2 /* 2 - e -(*+V2) 2 A 2 ] for | z | < L / 2 
Uz ~ [0 for \z\ > L/2 ' 

. (13) 

and u x = u y = 0. The variational equations give rise to 
a Schrodinger equation. The sign of A is nipped relative 
to the case of the contraction. In this case, A is every- 
where positive reflecting a repulsive effective potential. 
This difference in sign gives rise to an important differ- 
ence between expansion and contraction. In the case of 
expansion, the effective potential displays two peaks in- 
stead of two wells. In the absence of the two peaks, the 
problem reduces to that of a particle in an infinite poten- 
tial well model. The wavefunction for the unperturbed 
ground state is now given by 



L C ° S (Z Z 



(14) 



The energy of such a bound state in a box of size L is 



E 



2 



(15) 



Now, consider the perturbed state with the potential 
given by A = dV ■ u which reads 



A = 



2duo 



(z + ^ )e -(*+lr) 2 A 2 _(z_|) e -(*-#)7* 



L 



(16) 

We may approximate A near its maxima by delta func- 
tions. The maxima occur at z — ± We express 
A as 



A = duo 



(17) 



The first order approximation to the perturbed ground 
state energy trivially yields 



E' = E+ I r x H z = E+ 8 -^ S in k7: 



Ly/2' 



(18) 



Replicating the steps of Section III we find that the effec- 
tive transition temperature T^ 1 for the case of expansion 
is 



7r 2 c 8g?uo . kir 
— sin 



air 



W2 



In this the system is cooled down, the faces would 

become supersolid after the bulk crystal as we cool down 
the crystal. A plot is given in Fig. [3] for = 2 K, 
u = 0.5, 4^ = 0-02 K, and ^ sin = 1 K. As 
in our prior analysis of the contraction, the large value 
of uo in Fig. [4] is chosen to vividly illustrate the elastic 
distortion associated with an expansion. 

It is worth highlighting the origin of the difference be- 
tween the cases of edge contraction and expansion. Both 
cases have different divergences of the displacement field 
(and thus different local density profiles). The local mass 
or equivalently, the vacancy density is what couples to 
the supersolid order parameter. Note, in case of a su- 
perconductor it is the charge density that couples to the 
order parameter. Both the displacement field and the 
spatial gradient are odd under spatial reflection. In our 
case, V • u is even under spatial reflection (it reflects the 
scalar mass density) and the two cases are physically very 
different even though the spatial profile of the displace- 
ment fields in both cases are related by a minus sign (see 
Eqs.((5]) and @). 



V. DISLOCATIONS 

Below, we will present detailed numerical and vari- 
ational calculations of the local transition temperature 
due to an edge dislocation using the formalism that we 
have employed thus far in this work. For a discussion 
of dislocations in the quantum arena see, e.g., Ref. [24] 
An analysis analogous to ours was done by Tonei^ who 
reached similar conclusions as we have. Some time after 
we discussed this phenomenon^, Ref. [35] considered the 
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Figure 3: (Color online.) (a) The transition temperature along a slab is plotted in a contour map for different cases 
from top to bottom: a slab with contraction, a uniform slab, and a slab with expansion near its edges. Whenever 
the elastic deformations are present, the local supersolid transition temperature is altered by comparison to the 
uniform solid. Near the edges, where the elastic deformation is present, the supersolid transition temperature is 
altered: T c increases at the boundaries in the case of boundary contraction and decreases for an expansion near the 
boundaries. The dark (black) line is associated with contraction, red line with the normal crystal, and blue with the 

expansion. With the lattice constant to unity, the parameters k = 1, T~ = 2 K, uq = 0.5 and -j— 1 = 0.2 K. (b) We 
plot the effective transition temperature for an undeformed crystal (straight line), and that with a 
compression/expansion of its boundaries. 




Figure 4: (Color online.) The displacement field corresponding to an expansion near the edges. Plotted is the 
displacement field (black solid) given by Eq. ( 13 ) with L — 20, Uq — 0.5, d = 1 and k = 1. The effective potential (red 
dashed) is given by Eq. (16 1. The large (exaggerated) displacement highlights the expansion near the boundaries. 



problem of dislocation line filaments which become super- 
solid while the bulk is non-supersolid. This is markedly 
different from our approach where both the bulk and the 
dislocation core become supersolid at transition temper- 
atures that differ by small amounts. The small change 
in the ordering temperature is imperative in our pertur- 
bative approach of linearly expanding a(T) in Eq. ([!]) 



about the bulk supersolid transition temperature and in 
neglecting the cubic terms in Eq. ([3| when solving the 
effective Schrodinger type equation of Eq. Q . 

Many displacement fields can correspond to a given 
"Burgers vector" b describing a dislocation. The Burgers 
vector is defined by a circuit integral around a dislocation 
coreCSl = dg. Vu 7 ) for a large contour K around 



7 



the dislocation core) describing a dislocation. We will an- 
alyze one such particular set of displacement fields. All 
of these displacement fields are related to one another 
via a smooth deformation u — > u + v. Here, v is a non- 
singular vector field with a vanishing associated circula- 
tion: § c ds- Vv 7 = around any closed contour C. As we 
saw in the earlier sections, a smooth displacement field 
(corresponding to, e.g., a contraction or an expansion) 
can, on its own, raise or lower the effective local super- 
solid transition temperature. Thus, the effective change 
in T c , which we turn to next, will generally depend on 
the detailed form of the displacement fields u correspond- 
ing to a given dislocation. In what follows, we consider 
a particular minimal displacement field form that corre- 
sponds to symmetric unidirectional displacements about 
a lattice direction of an unstrained crystal. With x and y 
denoting the horizontal and vertical Cartesian directions, 
we consider a particular displacement field in Figs 



5(a)) 



and |5(b)j that corresponds to a dislocation with a Burg- 
ers vector b — be x . In what follows, the spatial extent of 
the dislocation core will be set by k. 

The following displacement field describes such an edge 
dislocation, 



t(x,y) = -— e {x +y >' sgn(x) cos 
2n 



(19) 

where sgn{x) is the sign of x, i.e., sgn(x) = [26 (x) — 1] 
with 8(x) the Heavyside function. In Eq. (19), the mag- 



nitude of the Burgers vector b cannot exceed the inter- 
atomic lattice spacing. Furthermore, realistically, k may 
be of order of 10 (lattice constants). We may derive an 
effective potential from the displacement in the same way 
we did for the above two cases (Eq. ([2])). In this case, an 
analytical solution to Schrodinger equation is not possible 
and we will resort to a numerical solution. The effective 
potential energy is provided in Fig. [6j 

We approximate the partial derivatives in the 
Schrodinger equation by finite differences and use the 
numerical Gauss-Seidel method for solving iteratively 
a system of linear equations in conjunction with over- 
relaxation. Our relaxation scheme shows that the wave- 
function localizes rapidly to the region near the origin 
where the dislocation core sits. Both the initial seed and 
the final numerical result are depicted in Fig. [7j The nu- 
merical solution to the Schrodinger equation illustrates 
that there is a change in transition temperature. 

We now turn to an approximate analytical solution. A 
contending variational state that is localized about the 
dislocation core is given by 



ip(x,y) 



-cr(x 2 +y 2 ) 



(20) 



From Eq. (19), we can compute the effective potential 



energy, 

V = X = gJ(V • u) 
bd 



2tt 



sgn(x)e 



x(x 2 + y 2 ) 



2x 



cos 



yjx 2 + y 2 



(21) 



The Hamiltonian H = [— cV 2 + V] corresponds to the 
Schrodinger equation of Eq. Q . The expectation value 



is 



(H) = 2ca - 



abd 

2Jnk 2 (2<T+±)3 



> E, 



ground ■ 



(22) 



An extremizing variational value of a is given by 



k 2 




(23) 



Substituting the above equation into Eq. (22), we obtain 

' Pk 2 - 1 ' 



(H) 



k 2 ' 



2/5 



bd 



4V7^ 3 fc 4 



(24) 



where /3 = ( 8 J^fc 4 ) ■ Similar to the earlier two 
cases of boundary deformations, the effective transition 
temperature can be found by approximating a(T) ss 
a(T — T®). The transition temperature 



Tf" = T? - -c 



bd 



\4^/3 3 k 4 



(25) 



In Fig. Hj we plot AT C = T^f - T c ° for different Burg- 
ers vectors b as a function of dislocation core size k for 
fixed values of the parameters in the GL functional of 
Eqs. and pi): c = d = 1. Depending on the choice of 
parameters an increase or decrease in the local transition 
temperature relative to a bulk transition is possible. As- 
suming a typical Burgers vector magnitude of b = 1 for 
an edge dislocation in an hep crystal, with parameters 
c = d = 1, a core size radius of R — k/y/2 < 3.5 lattice 
constants will result in a supersolid dislocation core prior 
to the bulk supersolid transition. 



VI. CONCLUSIONS 

In summary, we find that elastic deformations in super- 
solid lead to local changes in the transition temperature. 
For a positive coupling constant d in Eq. ([2| we obtain 
the results: 

1. Edge contraction increases the supersolid transition 
temperature at and near the edges. 

2. Edge expansion decreases the supersolid transition 
temperature at and near the edges. 



/ /r.\\\ 



Displacement of the Atom Planes 
Around the Origin 



(a)An edge dislocation 



(b)The displacement field corresponding to an edge dislocation 



Figure 5: A schematic of an edge dislocation. At left (a), shown are rows of atoms. The presence of an edge 
dislocation is manifest in the appearance of a different number of vertical rows of atoms above and below the 
terminal dislocation point. The corresponding displacement field is shown at right (b). 
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Figure 6: (Color online.) The effective potential energy 
of Eq. (21 1 corresponding to the displacement function 
of Eq. ( 19 1. In this figure, b = d = c = and k = 10. 
That is, /? = 0.0135 and a = 0.0017. The maximum 
strain is 0.086 (in units of the lattice constant). 



3. The local supersolid transition may be enhanced or 
suppressed near a dislocation core. 



This implies the observation of interesting effects. For 
example, for edge contractions, we would find that, be- 
low a certain temperature that is higher than the su- 
persolid transition temperature, a sample of supersolid 
would have its supersolid edges partially decouple from 
its bulk crystal. Of course, the effects of elastic deforma- 
tions on the supersolid transition are not limited to the 
few selected cases studied here. For example, the same 
physics applies to point defects like interstitials and va- 
cancies, as well as to extended defects like grain bound- 
aries and inclusions or voids. The above conclusions were 
based on the assumption of a positive coupling constant 
d in Eq. Formally, for negative d, our conclusions 
would have been inverted- an expansion would enhance 
the local supersolid transition while a contraction would 
reduce the supersolid transition temperature. 

Similar effects are found elsewhere in regions that lo- 
cally expand or contract. In Ref. [23] it was shown how a 
dislocation condensate may generally enhance and trig- 
ger superfluid behavior via a Higgs type mechanism. The 
presented GL approach of elastic deformations on the 
supersolid transition temperature is quite general and 
also applies to superconductors. In fact, dislocation de- 
fects and lattice-mismatched interfaces in superconduc- 
tors are known to create nonuniform strain and changes 
to the superconducting transition temperature, which 
have been studied extensively.^ Thus, our calculations of 
the changes in the local transition temperature due to a 
nonuniform elastic strain coupling in a Ginzburg-Landau 
approach are not limited to supersolidity and may as well 
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Figure 7: (Color online.) The supersolid order parameter associated with the effective potential of Fig. [6] Shown is 
the Gauss-Seidel solution of Eqs. ^ and Q with c = 1 for the dislocation profile of Eq. (19). Left: An initial seed 
state. Right: The final "wavefunction" (supersolid order parameter) ip. The localized bound state of the order 
parameter is evident. 



apply to superconductivity. 
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